We 
明寺 


7No.2 


i a ee A 
\ NaA IVD Ti 


pz 


第 37 卷 第 2 期 计算 机 应 用 研究 一 Vol 
录用 定稿 Application Research of Computers Accepted Paper 
第 四 方 物 流 协同 路 径 定 制 问题 研究 
符 妍 ', 黄 敏和 李 波 ) 
人 


发 


(1. 沈阳 工程 学 院 信息 学 院 , 沈阳 110136 2. 东北 大 学 信息 科学 与 工程 学 院 流程 工业 综合 自动 化 教育 部 重点 实验 
室 , 沈阳 110819) 


摘 要 : 为 了 解决 第 三 方 物流 (3PL) 供应 商 面临 的 高 运输 成 本 问题 ， 以 第 四 方 物流 (4PL) 供应 商 的 角度 ， 提 出 了 
多 第 三 方 物流 (3PL) 供应 商 协同 为 客户 定制 路 径 的 问题 。 针 对 该 问题 在 模型 求解 上 要 同时 考虑 路 径 与 3PL 供应 商 
的 特点 ， 设 计 了 基于 k- 短 路 的 混合 粒子 群 算法 (K-PSO)。 实 验 分 析 中 ， 通 过 开 -PSO 算法 与 遗传 算法 和 枚 举 算 法 对 不 
同 节点 数目 以 及 不 同 3PL 供应 商 个 数 下 仿真 算 例 的 计算 表明 了 算法 的 有 效 性 。 最 后 ， 通 过 对 3PL 转运 成 本 的 改变 ， 
展示 了 4PL 协同 运输 的 优势 。 
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Research on 4PL collaborative routing customization problem 
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Abstract: In order to solve the high transportation cost problem faced by the third party logistics (3PL) suppliers, from the 
fourth party logistics (4PL) supplier’s point of view, a multi 3PL Collaborative Routing Customization Problem is proposed. 
As the problem needs to consider the route and the 3PLs simultaneously, a hybrid Particle Swarm Optimization based 
K-shortest path algorithm (K-PSO) is designed. In the experimental analysis, the results show the effectiveness of the 
K-PSO by comparing with the GA and enumeration algorithm with different nodes’ and 3PL suppliers’ number. Finally, 
through the changes of the 3PL suppliers’ transit cost, it shows the advantages of the 4PL coordinated transportation. 

Key words: collaborative transportation; particle swarm optimization; fourth party logistics; k-shortest path algorithm 


0 引言 下 可 找 约束 最 短路 的 特征 ， 还 要 计算 转运 节点 的 转换 费用 和 
时 间 。 如 果 使 用 简单 图 求 出 所 有 路 径 ， 虽 然 理论 可 行 ， 但 解 

第 四 方 物流 (forth party logistics, 4PL) 协 同 路 径 定制 问题 空间 会 非常 庞大 。 因 此 ， 本 文 首先 利用 KK- 短路 算法 求解 简单 
(4PL collaborative routing customization problem, 4PLCRCP) 图 的 前 K 短路 径 , 然后 再 利用 粒子 群 算法 选择 路 径 上 的 供应 
是 指 单 第 三 方 物流 (Third Party Logistics，3PL) 供 应 商 很 难 胜 商 。 提 高 了 算法 效率 的 同时 增强 了 算法 的 有 效 性 。 


任 或 者 需要 付出 较 高 代价 才能 完成 客户 交 给 的 运输 工作 时 ， 鉴于 以 上 分 析 ， 本 文 基于 4PL 协同 运输 模式 ， 建 立 了 同 
4PL 供应 商 利用 其 拥有 的 信息 资源 协同 多 个 3PL 供应 商 共 ”时 考虑 3PL 供应 商 停靠 与 转换 成 本 的 4PLCRCP 数学 模型 。 
司 为 客户 定制 运输 方案 的 问题 。 设计 了 基于 K- 短 路 算法 的 混合 粒子 群 算法 (K-PSO)。 实 验 


有 关 协 同 物流 的 研究 ， 其 核心 任务 是 对 网 络 资源 进行 整 ” 分 析 中 ,针对 不 同 节点 数目 和 不 同 3PL 供应 商 数目 下 的 实例 
本 规 划 、 调 配 获取 和 和 集成 优化 由 。 对 于 4PLCRCP，4PL 供 应 。 计算 ， 以 及 与 遗传 算法 和 枚 举 算法 计算 结果 的 对 比分 析 表 明 
前 不 仅 要 协同 多 个 物流 分 支 ， 还 肩负 着 为 物流 分 支 选择 3PL 了 K-PSO 算法 的 有 效 性 。 最后, 通过 对 不 同 转 换 成 本 下 的 实 
供应 商 的 工作 。 在 之 前 的 研究 中 ， 任 意 两 节点 间 的 3P 例 计算 表明 :与 客户 只 使 用 单一 3PL 供应 商 的 运输 模式 相 比 ， 


Ld 
外 


有 有 被 认为 是 独立 的 中。 而 实际 运输 中 很 多 3PL 业务 范围 大 、 选择 4PL 协同 运输 模式 不 仅 可 以 增加 任务 按时 完工 的 比率 还 
只 是 在 个 别 路 段 存在 劣势 。 因 此 ，4PL 应 利用 其 拥有 的 信息 能 为 客户 节省 一 定 的 费用 。 
资源 为 客户 选择 3PL 供应 商 ， 从 而 使 配送 方案 更 具 优势 。 本 
可 是 描 壕 王 
在 4PLCRCP 中 ， 由 于 协同 机 制 的 加 入 ， 使 4PL 要 顾及 1 ”问题 描述 与 建 模 
各 转运 节点 的 具体 情况 。 许 多 学 者 对 4PL 问题 的 求解 方法 进 1.1 问题 描述 
行 了 研究 2-5。 但 是 在 4PLCRCP 中 ,问题 除了 在 多 重 图 假设 某 4PL 供应 商 拟 承 揽 一 项 由 供应 点 v 到 目的 点 的 
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运输 任务 ， 其 中 ~v 为 中 转 节 点 。 网 络 由 7 个 3PL 供应 商 
组 成 。4PL 可 以 根据 3PL 供应 商 的 运输 范围 指派 任务 。 每 个 
3PL 供应 商 在 运输 过 程 和 在 经 过 中 转 节 点 时 均 需 要 一 定 的 时 
间 和 费用 。 如 果 4PL 在 相 邻 路 径 使 用 不 同 的 3PL 供应 商 , 还 
会 产生 转换 3PL 的 时 间 和 费用 (起 始 和 终止 节点 视 为 转换 节 
磺 加 


将 上 述 问题 的 3PL 供应 商 看 做 边 并 以 编号 1,.…,r 标志 ， 
于 是 得 到 4PLCRCP 运输 网 络 多 重 图 GWV,E) 〈 其 中 ，Y 为 节 
点 集 ，E 为 边 集 )。 
基于 上 述 假 设 ， 为 建立 4PLCRCP 模型 需要 定义 的 参数 
及 变量 如 下 : 
1) 参数 说 明 
表示 运输 网 络 中 3PL 供应 商 的 数量 ; 
7 表示 交 货 时 间 ; 
Tx 表示 :3PL 供应 商 基 将 货物 由 点 i 运输 到 点 /的 时 间 ; 
7 表示 在 节点 i 的 停靠 时 间 ; 
7 表示 在 节点 ;更 换 3PL 供应 商 的 时 间 ; 
-1, i 是 始 发 节点 
b; :点 i 对 货物 的 需求 量 ,b=14 1， i 是 需求 节点 
0， 否则 
cx 表示 3PL 供应 商 基 将 货物 由 点 i 运 至 点 j 的 费用 ; 
C, 表示 在 中 转 节点 i 的 停靠 费用 ; 
GC, 表示 在 节点 i 更 换 3PL 供应 商 的 费用 。 
2) 变量 定义 
zx 表示 供应 商 k 是 否 承 担 i 到 j 的 运输 任务 , 即 
[ 供应 商 k 负 责 i 到 j 的 运输 任务 


IT 


次 二 10， 否则 

表示 节点 上 上 是否 承担 运输 任务 ， 即 
”_ 有， 点 i 承担 运输 任务 
“10， 否则 


5 表示 途径 点 i 时 ,是 否 在 点 i 更 换 3PL 供应 商 。 即 


上 人 y; 一 1 且 在 点 i 更换 3PL 供 应 商 


0， 否则 
其 中 ij ef{l...,n},ke{l,...,r} 。 
1.2 数学 模型 

本 文 目标 是 求 得 一 条 满足 约束 条 件 的 费用 最 短路 R 。 在 路 
径 R 下 ， 要 同时 考虑 3PL 供应 商 、 转 运 节 点 以 及 转换 3PL 供 
应 商 的 时 间 和 费用 。 根 据 上 述 描述 ， 建 立 的 模型 如 下 : 


Z=min( YD Cx t+ YOY + CC zi) (1) 
i=1 j=l k=l i=l i=1 
Sb 
DD Tx th Ty + Tz <T (2) 
i=1 j=l k=l i=l i=l 
DD a=b, jE{l2,...,n)} (3) 
i=1 k=l i=1 k=l 
DD = jE{2,..., n} (4) 
i=l k=l 
i iE{1,2,...,n-1) (5) 
j= k=l 
站 七 用， EE 位 2 二 (6) 
y=1 (7) 
y, =1 (8) 
x , is 5; = 0 或 1 ,b,j e{l,2,...,n)}, ke{l,2,...,7) (9) 


在 上 面 的 模型 中 ， 式 (1) 为 目标 函数 ， 表 示 在 运输 途中 、 
在 转运 节点 停靠 和 转换 的 费用 总 耗 ， 式 “2) 表示 到 达 目 的 地 
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的 时 间 不 大 于 客户 要 求 的 时 间 ; 式 (3) 表 示 流 量 平衡 ， 式 (4)(5) 
表示 选择 的 点 是 从 起 点 到 目的 节点 的 通路 ; 式 (6) 表 示 在 i 点 转 
换 则 一 定 在 i 点 中 转 ， 式 (7)(8) 分 别 表 示 路 径 的 起 始 和 终止 节 
点 ; 式 (9) 表 示 Nijt 1 和 Zi 为 0-1 变量 。 
2 ”基于 K- 短 路 算法 的 混合 PSO 算法 


2.1 主要 步骤 
在 K-PSO 中 , 求解 过 程 包含 初 始 解 生 成 和 解 的 更 新 两 个 
K 


阶段 。 在 第 一 阶段 ， 首 先 利用 K- 短 路 算法 求 前 K 短路 径 ， 
然后 根据 种 群 规模 等 比例 扩展 K 短路 数目 , 最 后 随机 生成 路 
径 上 的 3PL 并 计算 适 值 。 在 第 二 阶段 ， 使 用 第 一 阶段 的 结果 
作为 初始 种 群 ， 然 后 利用 PSO 算法 选择 路 径 上 的 3PL。 每 代 
种 群 中 ， 对 上 代 KK 短路 的 平均 适 值 进行 从 小 到 大 排序 。 然 后 
利用 旋转 轮 盘 法 中 确定 第 K 短路 的 数目 。 基 于 上 面 描 述 ， 算 
法 的 主要 步骤 如 下 : 
a) 利 用 K 短路 算法 求 前 K 短路 径 。 根 据 KK 短路 的 结果 确 
定 开 值 。 等 比例 扩展 种 群 规模 后 ,初始 化 路 径 上 的 3PL 供应 
商 并 计算 适 值 ; 

b) 根 据 K 短路 的 平均 适 值 由 轮 盘 法 分 配 第 KK 短路 径 在 种 
群 中 的 个 体 数 目 ; 
c) 使 用 PSO 算法 求 路 径 上 的 3PL 供应 商 并 计算 适 值 。 
d) 对 于 不 同 的 K 短路， 是否 有 更 好 解 ? 是 : 更 新 该 KK 短 
路 的 最 好 解 ， 否 则 ， 不 变 ; 

e) 是 否 达 到 最 大 迭代 次 数 ? 是 : 转 Step6; 否 : 转 Step2; 

f) 算 法 运行 结束 ， 比 较 各 个 K 短路 的 最 好 解 ， 输 出 最 优 


者 。 
2.2 编码 结构 及 初始 解 

根据 问题 特点 ， 解 的 编码 结构 由 路 径 上 的 节点 外 
PATH =[p,…,p1] 和 边 集 PL=[pl,…,pli1] 两 个 部 分 组 成 。 其 品 
Pp;( pieV ) 表 示 路 径 上 的 第 i 个 节点 ， pl (pl eB) 表示 使 用 上 


第 i 个 3PL 供应 商 (ieE{l,...,n} )。 


在 初始 解 生成 阶段 ， 路 径 上 的 节点 集 PATH 由 K 短路 算 
法 得 到 。 等 比例 扩展 节点 集 的 数目 至 种 群 规 模 ， 通 过 随机 生 
成 3PL 供应 商 的 方法 构成 解 的 边 集 , 从 而 构成 解 的 初始 种 群 。 
2.3 解 的 更 新 
在 解 的 更 新 阶段 , 首先 利用 PSO 算法 的 速度 位 置 更 新 式 
中 计算 路 径 上 3PL 供应 商 的 取 值 (0 或 者 1)， 然 后 使 用 轮 盘 
赌 方法 确定 3PL。 具 体 方案 如 下 : 
假设 第 8 个 粒子 位 置 : PL = Cnxp， …,Xow) ， 速 度 : 


DD Ar 


eT 


Vy = Von Vir)， 其 中 8=1…PS ，PS 表示 种 群 规模 ， 
7 n} 。 
VkE{L2,.…,7), 第 8 个 粒子 的 速度 Vox 和 位 置 x 更 新 公 
式 如 下 : 
Va Vi + CE(PS, i) + CMP- Xi ) (10) 
un | random < S(Ve) GD 
【0, 其 它 
4 [1 轮 盘 赌 选中 
多 (2) 
0 其 他 
其 中 : SVA)=1/(+exp(-V4)) ，d =1,2,.….,GEN ， GEN 表示 迭代 


次 数 ， Gl 和 C> 为 学 习 大 于， 和 7 为 [0,1] 的 随机 数 ， ps 为 粒 
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子 8 所 经 历 的 最 好 点 ， Ps 为 种 群 内 所 有 粒子 经 历 的 当前 最 
好 点 。 
根据 式 (10) 产生 的 V (k=142,.….,r )， 由 式 (11) 生 


i 到 j 的 3PL 的 0-1 取 值 。 分 别 计 算 3PL 取 值 为 1 时 的 路 

适 值 ,然后 按照 从 小 到 大 的 顺序 进行 排序 。 设 参数 4Es(0.D ， 

导 到 基于 序 的 评价 函数 : 

eval(k) = 1(1- 1)" (13) 

选择 过 程 以 旋转 轮 盘 的 方法 为 基础 ,根据 评价 函数 (13)， 
旋转 轮 盘 并 由 式 (12) 确定 最 终 被 选中 的 3PL 供应 商 。 

2.4 适 值 函数 


奏 芯 逮 
二 


利用 K- 短 路 算法 求解 K 短路 径 时 ， 未 考虑 路 径 的 时 间 
约束 。 本 文 将 时 间 约 束 作为 罚 值 加 入 到 目标 函数 中 ， 从 而 形 


成 适 值 函数 式 (14)。 
| TT>T 
[Zz 其 他 


其 中 g(a 三 DD) 为 惩罚 系数 。 
3 ”数值 实验 


3.1 实验 数据 

针对 4PLCRCP 的 特点 ， 假 设 网 络 节点 集 V 中 的 n 个 节 
点 在 100*100 范围 内 随机 生成 -NO 表示 可 能 与 节点 i( vieV ) 
有 边 相 连 的 邻接 点 个 数 。D 为 邻接 矩阵 (其 中 DGi, j)=r 表示 
i 到 j] 有 r 条 边 ，i,jeV )。d(i, j).costtiD 和 d(i, j).time( 让 分 别 
表示 节点 i 到 节点 j 的 第 六 (ii=1,...D 条 边 的 费用 和 时 间 属 
性 。 实 验 数据 分 别 由 下 面 三 个 阶段 生成 : 

a) 生 成 与 节点 i 相连 的 节点 。 

VieVUjznUjzi， 如 果 |j-i|<NO， 随 机 产生 1~NO 之 间 
的 随机 整数 世 。 如 果 t>NO-1, 则 设 定 i 到 j 之 间 没 有 边 相 连 ， 
即 DGi,j)=0; 否则 假设 i 到] 有 r 条 边 相 连 ， 即 DGi, j)=r。 

b) 生 成 边 (3PL 供应 商 ) 的 费用 和 时 间 属 性 。 

Vi,jeV ，VYiier， 如 果 D(i,j) *0, 随机 产生 0~1 之 间 的 
随机 数 世 。 假 设 当 t2>0.2 时 ,i 到 j 第 让 条 边 的 运输 费用 为 
1~15 之 间 的 随机 整数 ， 即 d(i,j).cost(ii)=rand(1,15); i 到 j 第 
让 条 边 的 运输 时 间 为 两 点 之 间 的 距离 除 以 费用 ， 即 
dGj).time(iD=|-jlydGjj).costtiD。 当 t2< 0.2 或 对 于 其 他 节点 ， 
3PL 供应 商 的 费用 和 时 间 为 无 穷 大 。 

c) 生 成 节点 的 费用 和 时 间 属 性 。 

vieV ,对 节点 i 的 中 转 和 转运 费用 分 别 赋 予 2~12 的 随机 


整数 ;对 节点 i 的 中 转 和 转运 时 间 分 别 赋予 4~14 的 随机 整数 。 


3.2 结果 分 析 
腿 设 某 4PL 供应 商 拟 承 揽 一 项 


(14) 


节点 =1 到 目的 点 


%=8，7=5 的 运输 任务 。 数 据 由 3.1 节 NO=3 时 产生 。 算 法 
使 用 Matlab7.0 并 在 core2 2.83GHz Pc 上 实现 。 经 多 次 实验 ， 


较为 合理 的 算法 参数 设置 如 下 : 

K=3,PS=5,GEN=20,c1=c2=2, § =0.5, 7 =0.8, 4 =0.5, w=100 

段 设 T=26 天 ， 得 到 的 最 优 解 为 [1, 3, 6, 8] [1, 5, 5] 

这 意味 着 ， 当 决策 者 (4PL) 需要 在 26 天 内 将 货物 由 节 

点 1 运 至 节点 8 时 ， 使 用 编号 为 1,5,5 的 3PL 供应 商 分 别 负 
节点 1 到 3、 节 点 3 到 6 和 节点 6 到 8 的 运输 工作 。 这 一 

运 桨 过 生 站 省亲 25 天 ， 费 用 157.7。 

图 1 展示 了 n=8 时 K-PSO 算法 的 收敛 情况 。 

为 了 增强 算法 求解 4PLCRCP 的 说 服 力 ， 随 机 产生 
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过 改进 文献 [2] 设 计 的 遗传 算法 (简称 GA) 和 枚 举 算法 《〈 简 
称 EM)， 分 别 对 上 述 算 例 进行 求解 。 假 设 每 个 算 例 运行 50 
次 ， 得 到 的 结果 如 表 1 所 示 。 


280 
二 
*# 解 的 均值 

260| ， @ 当前 最 好 解 
最 240| “ ] 
好 
解 220| +** 
收 
伍 200| ” 
情 
、 站 六 本 
况 180 Be 

ooo。*++ 
160 上 Co oo 四 由 由 由 四 帅 电 帅 电 遇 遇 遇 晶 遇 各 各 各 
140 
10 15 20 25 30 
迁 代 次 数 
图 1 8 节点 K-PSO 算法 收敛 曲线 
Fig. 1 Convergence curve of the 8-node K-PSO algorithm 
表 1 中 ,“-” 表 示 在 相应 算法 中 不 存在 该 项 信息 。 对 于 


EM， 当 n=8 时 ， 总 计算 时 间 79 秒 。 这 一 结果 与 使 用 K-PSO 
和 GA 计算 得 到 的 结果 一 致 。 当 n=30 时 ，EM 无 法 在 有 限时 
间 内 求 得 最 优 解 。 而 K-PSO 算法 在 最 好 解 的 均值 、 偏 差 等 方 
面 均 优 于 GA。 男 一 方面 ，GA 的 种 群 规模 远大 于 K-PSO 的 
种 群 规 模 。 造 成 这 种 现象 的 原因 是 协同 机 制 的 加 入 使 得 GA 
很 难 顾及 每 一 转运 点 , 从 而 计算 的 稳定 性 要 依赖 于 种 群 规模 。 
当 n=50 和 n=100 时 ， 两 种 算法 在 偏差 、 最 坏 解 和 均值 等 方 
而 的 差异 更 加 明显 。 因 此 ， 从 表 1 可 以 看 出 ,虽然 K-PSO 因 
计算 K 短路 耗费 了 一 定 的 时 间 ， 但 K-PSO 算法 求解 的 平均 
质量 要 优 于 GA。 

表 1 K-PSO 与 GA 和 EM 的 对 比 情况 


Table 1 The comparison of the K-PSO with GA and EM 
algorithm V E 7T K PS GEN Best Bad 418 msd Time (s) 
K-PSO 8 85 26 3 5 20 157.7 157.7 157.7 0 <1 
GA 8 85 26 - 20 0 15377 137,7 13577 和 <1 
EM 8 85 26 - - - 157:7 - - - a 


K-PSO 30 ‘B13. 在 起 10 50 249.6 259.6 250.8 2.5 3 
GA 30 815 110 - 70 40 249.6 261.4 254.6 5.7 1 
EM 30 815 110 - 
K-PSO 0 1420 1553 3 1s 70 349.6 361.6 351.5 3.8 8 
GA 50 1420 155 - 130 60 349.6 378.5 355.4 11.2 3 
EM 50 1420 155 - 
K-PSO 100 3160 280 8 25 100 587.1 611.2 590.4 5.2 14 
GA 100 3160 280 - 220 100 587.1 635.6 596.9 17.9 污 
EM 100 3 160 280 - - - 本 e 本 

图 2 展示 了 n=30，NO=5，r=1, 5, 10, 15, 20 时 ，GA 和 
K-PSO 计算 偏差 的 对 比 情况 。 
当 r=1 时 ， 任 意 两 点 之 间 只 有 一 个 3PL， 于 是 4PLCRCP 
转换 为 带 有 约束 的 最 短路 问题 。 随 着 3PL 数目 增多 ， 两 种 算 
法 的 偏差 有 所 增加 。 但 是 当 r=10, 15 和 20 时 ，K-PSO 的 求 
解 偏差 缺 有 一 定 幅 度 的 减少 。 产 生 这 种 现象 的 原因 是 随 着 
3PL 个 数 增多 , K-PSO 算法 的 选择 余地 变 大 ,而 GA 却 因 3PL 
数目 增加 导致 算法 求解 能 力 变 弱 。 
3.3 问题 分 析 

在 4PLCRCP 中 ，4PL 的 决策 结果 不 仅 与 3PL 供应 商 的 


NO=5、n=30，NO=8、n=50 以 及 NO=8、n=100 的 算 例 。 通 


运输 时 间 和 费用 有 关 ， 还 与 节点 的 转换 成 本 有 着 干 丝 万 缕 的 
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联系 。 表 2 展示 了 n=30 时 4PL 协同 模式 与 单 3PL 供应 商 运 
输 时 的 对 比 情况 。 在 协同 模式 下 ， 当 C =T =0 时 ，3PL 供应 
的 转换 成 本 为 0， 这 时 4PLCRCP 转换 为 4PLRP 问题 外 。 
随 着 转换 费用 增加 ，4PL 尽量 使 用 相同 3PL 供应 商 。 与 此 同 
对 ， 随 着 转运 时 间 增 加 ， 运 输 方案 也 会 随 之 改变 。 从 表 2 可 
以 看 出 ， 在 满足 时 间 约 束 的 条 件 下 ， 相 同 转换 费用 不 同 转 换 
对 间 下 的 运输 费用 相同 。 但 是 随 着 转换 时 间 的 增加 ， 如 
C=1,7=20 ，4PL 便 无 法 按时 完成 任务 ， 这 时 总 的 运输 费用 
也 会 明显 增 大 。 由 于 3PL 业务 范围 的 限制 ， 在 n=30 的 算 例 
中 ， 只 有 3PL 供应 商 1,2,5 能 够 独自 完成 部 分 任务 。 如 果 使 
用 它们 单独 运输 ， 虽 然 不 需要 支付 转换 成 本 ,但 可 能 会 增加 
总 的 运输 费用 。 对 表 2 中 不 同 转换 成 本 下 完成 任务 的 个 数 ( 图 
3) 以 及 与 协同 运输 相 比 运输 费用 的 增加 幅度 (图 4) 上 来 看 ， 
协同 运输 模式 不 仅 可 以 节省 费用 ,还 在 完成 任务 个 数 〈 图 3) 
以 及 完成 任务 时 的 费用 上 (图 4) 明显 优 于 客户 使 用 单 3PL 
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五 


五 
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供应 商 的 结果 。 
20 -一 


02 4 6 8 10 12 14 16 18 20 
3PL 个 数 
图 2 不 同 3PL 数量 下 的 偏差 对 比 情 况 
Fig.2 The deviation under the different 3PL quantities 


表 2 不 同 运输 方式 下 转换 成 本 对 最 优 解 的 影响 


Table 2 The influence of the conversion cost on optimal solution under different transportation modes 


. 协同 运输 3PL(1) 3PL(2) 3PL(5) 
| | 路 径 3PL 供应 商 时 间 / 费 用 时 间 / 费 用 时 间 / 费 用 时 间 / 费 用 
0 0 [1511 18 22 23 30] [532323] 66/227 95/227 86/234 88/296 
1 [1 5 11 18 22 23 30] [111111] 95/229 95/229 84/236 88/297 
1 5 [1 5 11 18 22 23 30] [111111] 105/229 105/229 94/236 98/297 
1 10 [1 5 11 18 22 23 30] [222222] 104/236 115/749 104/236 108/297 
1 15 [1 5 9 16 23 30] [55533] 105/270 125/1 749 114/696 118/1 127 
1 20 [1 5 11 18 22 23 30] [222222] 124/1 696 135/2 749 124/1 696 128/2 127 
2 0 [1 5 11 18 22 23 30] [111111] 95/231 95/231 84/238 88/299 
2 5 [1 5 11 18 22 23 30] [111111] 105/231 105/231 94/238 98/299 
2 10 [1 5 11 18 22 23 30] [222222] 104/238 115/751 104/238 108/299 
2 15 [1 5 9 16 23 30] [55533] 105/273 125/1 751 114/698 118/1 129 
2 20 [1 5 11 18 22 23 30] [222222] 124/1 698 135/2 751 124/1 698 128/2 129 
3 0 [1 5 11 18 22 23 30] [et 95/233 95/233 84/240 88/301 
3 5 [1 5 11 18 22 23 30] [111111] 105/233 105/233 94/240 98/301 
3 10 [1 5 11 18 22 23 30] [222222] 104/240 115/753 104/240 108/301 
3 15 [1 5 9 16 23 30] [55533] 105/276 125/1 753 114/700 118/1 131 
3 20 [1 5 11 18 22 23 30] [222222] 124/1 700 135/2 753 124/1 700 128/2 131 
4 0 [1 5 11 18 22 23 30] [111111] 95/235 95/235 84/242 88/303 
4 5 [1 5 11 18 22 23 30] [111111] 105/235 105/235 94/242 98/303 
4 10 [1 5 9 16 23 30] [555514] 104/242 115/755 104/242 108/303 
4 15 [1 5 9 16 23 30] [55533] 105/279 125/1 755 114/702 118/1 133 
4 20 [1 5 11 18 22 23 30] [222222] 124/1 702 135/2 755 124/1 702 128/2 133 
18 
16 间 1600 
14 :| 一 1400 
供 
完 12 应 1200 - 
10 1000 
数 “ 
4 ] EF 中 
0 T 4 2 3 4 5 6 7 8 9 而 本 et 
协同 运输 [=1 [=2 [=5 转运 成 本 〔 费 用 ,时 间 ) 


图 3 不 同 运输 方式 下 4PL 完成 任务 的 情况 
Fig.3 The 4pls’ task under the transportation modes 


图 4 不 同 转运 成 本 下 3PL 供应 商 1,2,5 的 费用 增加 幅度 图 
Fig.4 The different transfer costs’ chart of the 1,2,5 3PL supplier 
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第 四 方 物流 协同 路 径 定制 问题 (4PLCRCP) 是 一 类 复杂 
的 优化 问题 。 本文 在 考虑 转运 成 本 基础 上 , 建立 了 4PLCRCP 
的 数学 模型 并 设计 了 基于 K 短路 的 粒子 群 算法 (K-PSO)。 
通过 K-PSO 与 GA 和 枚 举 算 法 的 比较 验证 了 K-PSO 的 有 效 
性 。 另 外 ， 实 验 结果 表明 : 如 果 4PL 能 够 有 效 整合 多 个 3PL 
供应 商 ， 那 么 即使 在 中 转 节 点 收取 转换 费用 的 情况 下 ， 仍 然 
可 以 起 到 整合 优势 3PL、 节 省 费用 的 目的 。 
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